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The calculation of the two-loop contribution to the Lamb shift in hydrogen is one of 
the most challenging projects in bound state QED IHJ^. Since direct numerical calculations 
with the use of Dirac-Coulomb propagators have not yet been completed, one has to rely on 
the Za expansion: 



= m ' (Z 1^40 + {Z a) S50 



ln='(Z a) -2 + ln'(Z a)"' B^2 + ln(Z a)"' Egi + B, 



60 



+... . (1) 



The leading order correction B^^ can be obtained from the slope of the electron form-factors 
Fx and -F2 at = 0. It is known analytically and its numerical value is quite small (for 
S-states including vacuum polarization) 

= 0.538941 . (2) 

The calculation of the next order correction B^q was completed only a few years ago inde- 
pendently by two groups in and @]. The value was surprisingly large 

= -21.5561(31) . (3) 

Moreover, this correction led to strong disagreement in He"*" Lamb shift with the most 
precise experimental value in while for hydrogen Lamb shift, led to agreement with 
the Mainz value for the proton charge radius 0. This large value of i?5o compared to i?4o 
indicates a very slow convergence or even might suggest a nonperturbative behavior of the 
two-loop contribution. Indeed, the direct numerical calculations of one diagram, the loop- 
by-loop electron self-energy by Mallampalli and Sapirstein in 0, shows that the value of 
this correction at Z = 1 is of different sign and magnitude, that the one based on the first 
two terms of analytic expansion. Moreover, this numerical calculation was in disagreement 
with the analytical value of -Bes in while it was argued in that this correction 
comes only from this one diagram in the covariant gauge. A year later another group 
calculated numerically this one diagram and found agreement with the analytic expansion 
including ln^(Za)~^ term. While this may suggest that the first numerical calculations 
were incorrect, a very recent, third numerical result by Yerokhin in [|l^ confirmed the first 
one [0. So, this situation with the two- loop contribution is very unclear. Moreover, the 
analytic value ln^(Za;)~^ term corresponding to all diagrams, was confirmed independently 
by several groups, so this situation is even more confusing. It was argued, by Yerokhin 
in |T0|, that the ln^(ZQ;)~^ term for this one loop-by-loop diagram is different from the 
total value of Sea, and in fact found an additional contribution. However, the value for 
this term, coming from all diagrams might be correct, because other diagrams may contain 
compensating terms. The goal of this work is to shed some light into higher order two-loop 
corrections and calculate all logarithmic terms: Bq3,Bq2, and -Bgi- We find that indeed the 
two-loop contribution has a very peculiar behavior, as the higher order term Bqi dominates 
and reverses the sign for the overall logarithmic contribution. In next sections we present 
some details of this calculation. First, a simple example is worked out to demonstrate the 
method, then we pass to the most difficult two-photon-loop diagrams and complete with 
remaining diagrams containing an electron loop. Conclusions with prospects of calculation 
of Bqq summarize this work. 
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I. SIMPLE EXAMPLE 



The example to demonstrate the calculational method is the asymptotic expansion of 

1 



P 



E-{H + uj[ 



■PI 



1^2 4^2 ^ 4- 12 ln(2) +4 ln(u;) 



9 5' Q 
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(4) 



(5) 



around large lo for ground state of the hydrogen atom. More precisely, we concentrate on 
the a;~^ term. For simplicity, we put here m = 1, a = 1. From one side -P(a;) is known 
analytically |Tl[ 

2Fi(4,2-r,3-r,C), (6) 



:i + r)8(2-r) 



where 
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so one could get this coefficient from here. However, our final goal is to calculate the two- 
loop contribution, for which no analytic formula has been derived so far. Therefore, we use 
a different approach, based on the effective Hamiltonian. First, we regularize the Coulomb 
interaction by the following replacement: 

\/(r) = -i^-l(l-e-^0- (8) 
With the regularized potential [P Pr) one can expand Pr in (H — E)/uj which leads to 

1 , , . 1 



Pr = - 



\p{h-eYp\ 



\v'{ry 



2 A + 8 ln(3) - 8 ln(A) - 2 , 



(9) 
(10) 



where the last expectation value is taken from ||T2[. The remaining part, which was left out 
by this replacement, is obtained from the subtracted forward scattering amplitude. Two 
photon exchange is 



P2 



(2^ 
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p p 
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'4, (11) 



where we keep only the uj~^ term {uj~^ and uj~'^ are subtracted out before the integration). 
The three photon exchange requires more subtractions. One Coulomb exchange between 
photon vertices gives P^a 
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Coulomb exchanges out of photon vertices gives P^b 
Ar = -2048 TT^ i F I F I 



(2 7r)3 J (2 Tif \p'^ p'^ + 2uj 
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2-32 ln(2) + 16 ln(3) 



(13) 



There is an imphcit subtraction at p' = for removal of small p' divergence. It corresponds 
to subtraction of lower order contributions. Additionally, only the uj~^ term is selected. The 
sum 

4-12 ln(2) + 4 ln(w) , , 

is independent of A in the limit of large A and agrees with that from the expansion of analytic 
formula in Eq. (^). The advantage of this method is the direct application to the two- loop 
Lamb shift. 



II. TWO-LOOP LAMB SHIFT 

The calculations of the two-loop Lamb shift in the order of a'^{Z is more complicated 
due to the presence of powers of ln{Z a). It reflects the fact that several energy and momen- 
tum regions contribute. For these calculations we introduce a number of cutoff parameters 
to separate different regions and calculate them independently. In Fig. |l| the integration 
region of two photon energies ui and U2 is split with the help of ei, e2, e'^, €3. Additionally 
A 'splits' the integration over electron momenta. The splitting itself, does not help too 
much. The key trick is the assumption that after expansion in Za one goes to the limits 
62 0,ei —>■ 0, €2 — * 0, e[ — > 0, A ^ 00, in the order as written. The two-loop contribution 
is split accordingly 

AE = El + Em + Ef + Eh, (15) 

and calculated separately, each term in the most convenient gauge. In the following sec- 
tions we calculate all logs. The constant term Bqq is left unevaluated, however we lay the 
groundwork for its calculation. 



III. CONTRIBUTION El 

Diagrams in the Coulomb gauge in NRQED are presented in Fig. ^. We calculate them 
first, for photon energies inside a rectangular box < ti;i < ei,0 < a;2 < 62,62 << ei, 
and after combine to the region Ei as shown in Fig. |l|. The expression derived from 
nonrelativistic QED for all these diagrams is: 
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(16) 



It is a two-loop analog of Bethe logs. We have not found a way to calculate its matrix 
elements analytically in a compact form, therefore we proceed in a different way. One finds, 
that Sl as in Eq. ([T6|) depends on a only through ei and €2- 
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£1 
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To find the logarithmic dependence, we differentiate over ei and €2 which with the help 
of €2 << ei leads to a much simpler expression. The first derivative leads to 



2a 



2 «2 .1 

duJ2 UJ2 S^s^r) (0 b* 77}- ^ 



(18) 



where 5w53(r) denotes first order corrections to 0, H, E due to 7r5^(r) operator. This integral 
was considered and calculated in the context of hyperfine splitting in hydrogen-like systems 
TB| , since the Fermi spin-spin interaction is also proportional to S^{r). The result from that 



paper which is extended here to any value of principal quantum number is: 



2a 
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E-{H + uj] 



■P 



a 9 F(n) 
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F{n) = — ln^e + Ine 
3 



2(l-2 1n(2)) + |(| + -i^-i-lnH + vI/H + C 
6 \A An"^ n 



where N has been calculated only for n = 1 

N = N{1) = 17.8299093, 



N{n), 
(20) 

(21) 
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and \1/ = r'/r with Euler F function and Euler C constant 



+ 



n — \ 



-C 



(22) 



We have introduced here a notation e = e/a^, which is to be used throughout this work. 
The result for n = 1 with S = m{a/'n')'^ is: 



In2(e2) + 2 (1 - 2 In 2) ln(e2) + N 



dei 3 

The second derivative, over €2, is httle more difficuh to calculate: 



(23) 



€2 



deo 



2a 
3 TT 
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One splits it into two parts, with the assumption €[ << €2- The first term A has the same 
form as that in Eq. ( p3D with €2 replaced by €[. The second term B is in turn split into two 
parts B = Bl + Bh, where Bl is calculated with the regularized Coulomb potential, as in 
Eq. (W). One can expand here in the ratio {H — E)/uj which leads to the expression: 



ei 



(25) 



energy levels in 



Both terms in above braces have already been calculated in the context of positronium 



2 3 



4 2 



2 

n 



+2(ln(n) - ^(n) - C) 



-4 + A-4 + 6 1n? 



(26) 
(27) 



with n = 1 in our case. Bh is the difference between B and -B^. In this difference only large 
electron momenta contribute, therefore it could be obtained in the scattering amplitude 
approximation, in the same way as P2 and P3 in a simple example in the previous section. 
The result is 



B 



H 



£ 



5 tt" - In f ^"j + 2 A In f ^"j - 50 ln(2) In (^) + 18 ln(3) In (^j 



+ln 
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(28) 



The complete B term is 



B = ^-[8 + 5 7r2 + 3 In(ei) - 6 ln(2) In(ei) - 2 In(ei) ln(e2) + ln(e2)^ - 3 ln(e;^ 
9 



+6 ln(2) ln(e;) + ln(e; 



(29) 
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We can now go back to Eq. (^4]) for the second derivative of £l which is a sum of A and B 



€■2 



^-[8 + — + 5 7r^ + 3 In(ei) - 6 ln(2) In(ei) - 2 In(ei) ln(e2) + ln(e2 



- n2i 



de2 9 ' 2 
The expression for £i which matches both derivatives is: 



(30) 



^L(ei,e2) = S 
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8 ln(2) In(ei) ln(e2) 4 In(ei) ln(e2)^ ^ 41n(e2) 
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(31) 



The constant term (no logs) is not included here. El as shown in Fig. ^ is integrated over 
the region which is a combination of three rectangles: 



El — £l{ -1^, ^] + £l \ -^i 



£j 



(32) 



IV. CONTRIBUTION Em 

In the one-loop case, contribution to energy, coming from photon energies > e is 

5E= (01^10), (33) 



Em is a correction to the Bethe log: 

2a 



E 



M 



6v(e^) / duju{(j)\p' 



?>'K^'''^" Jo ~" " E-{H + uj, 
It has the same form as Eq. (p3|), so after symmetrization ei ^ £2 it is: 
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V. CONTRIBUTION Ep 

Ep is the two-loop contribution with regularized Coulomb interaction and with both 
photon energies limited from below by e. It is a sum of three terms 



Ep = Ep + Ep + Ep , 



(37) 



defined and calculated as follows. Ep is a second order correction coming from V{ei) and 
V{e2) with V defined in (P^), here additionally with A-regularization 



E'p = m{e^] 



{E-H) 



(38) 
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The corresponding matrix element is given in Eq. (^), so Ep becomes 

s (it - 5 '"'^ (y - 5 '"P-O (-4 A - 16 m ^ - 4 - 64 14) , m 

One needs only In A term, since others do not give In a. Ep is the contribution from electron 
formfactors and E2 at = on relativistic (Dirac) wave function. We know it from 
the one-loop case that for vacuum-polarization Aqi = A40/2. The same holds for two-loop 
contribution, thus we have 

El = S\na-'^. (40) 

Diagrams with closed fermion loop are automatically included in the above formula. Other 
contributions coming from these diagrams are calculated in Section VII. 

Ep is the contribution from F[' and F2 calculated with nonrelativistic wave functions. It 
leads to the matrix element {(pl'V'^ 6^{r)x\(j)) which does not lead to In A. Hence, it does not 
contribute to In a. 



VI. CONTRIBUTION Eh 

Eh is the contribution obtained from the two-loop three-photon exchange forward scat- 
tering amplitude. It requires subtractions of terms, contributing to Lamb shift at lower 
orders. After subtractions it is finite and depends on ei,e2 and A = Aa. When combined 
with El and Ep, the dependence on €1,62 and A should cancel out. Having this in mind, 
the In a contribution could be obtained by the replacement A 1/a in Ep in Eq. (^9|). 
However, the constant term Bqq requires complete calculation of Eh, which we think is the 
most difficult of the contributions. 



VII. DIAGRAMS WITH CLOSED FERMION LOOP 

There is a small logarithmic contribution coming from diagrams with a closed fermion 
loop. They are partially included in Ep. Two other contributions Eyp, Eyp are the follow- 
ing. The second order correction coming from the one-loop vacuum polarization is 

E'vp = ^{-^y ('/'I^aW (^^^ SlirM £: (^)' In « . (41) 

The second contribution Eyp is electron self-energy in the Coulomb potential including 
vacuum polarization correction. It is calculated in the similar way, as previous corrections. 
One splits it into three parts 

EI.p = Cl + Cm + Ch. (42) 
Cl is a v.p. correction V = —(4/15) 6^{r) to the Bethe log: 
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Cm is a second order correction coming from self-energy and v. p. 
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(43) 
(44) 

(45) 
(46) 



Ch is given by the scattering amplitude. Since we calculate only the logarithmic part, instead 
of calculating Bh we replaced In A by — In a in the equation above. The logarithmic part of 
electron self-energy in the Coulomb potential including vacuum polarization correction is 



15 



- ln^a"2 + 4( g +ln2) In a 



(47) 



This completes the treatment of two-loop logarithmic correction 



VIII. SUMMARY 



The sum of all logarithmic terms in Eqs. ( ^^^pj^Tl) is 
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27 




Bq2 = 
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10800 
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= -0.296296, 

= -0.461891, 

55 TT^ 616 In 2 



27 135 
= 50.309654. 
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First of all the result for Bqi is surprisingly large, and reverses the sign of the overall 
logarithmic contribution. i?63 agrees with the result obtained first in p]. However, as it was 
pointed out by Yerokhin |1^, the loop-by-loop diagram is the source of additional terms, 
which were not accounted for in the calculation in p|. An additional result of this work 
is the state dependence of B coefficients which is obtained from n-dependence of matrix 
elements in Eqs. (p0|,[26|j2l 
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56i(n) 
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4 /304 
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n J 



(51) 



f 1^(2)) + ^ - - - ln(n) + nn) + c) . (52) 



n 



n-dependence of 562 agrees with the former result in |[T^ (apart from the misprint in the 
overall sign there). -Bgi depends on A^-coefficient, the Dirac delta correction to Bethe logs. 
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which has not been calculated yet for other states than IS, therefore its complete state de- 
pendence is unknown. However, one may expect to a good approximation is independent 
of n, as it is for Bethe logs. 

Because of the large value of Bqi theoretical predictions for hydrogen Lamb shift are going 
to be changed. The total logarithmic contribution is 16.9 kHz for the IS state, compared 
to the previous one, based only on Bq^ -28.4 kHz. Theoretical predictions for Lamb shift in 
hydrogen with proton radius Vp = 0.862(12) fm from W^, using recent updates: analytical 



calculations of the three-loop contribution by Melnikov and Ritbergen in |]T6[ and direct 
numerical calculation of one-loop self-energy by Jentschura et al. in [^] are (see details in 
the appendix) 

Ei(15)th = 8 172 816(10) (32) kHz, (53) 
El{2S-2P„2U = 1057 842(1) (4) kHz, (54) 

where we assumed for B^q = 0±100, which gives the first uncertainty. For P-states we neglect 
S-terms completely. The second uncertainty comes from the proton charge radius. Since 
it dominates the theoretical error, we emphasize the importance of the muonic-hydrogen 
measurement, from which could be precisely obtained. Current theoretical predictions 
agrees well with the most precise experimental values: 

^l(15)cxp = 8 172 837(22) kHz [18, m, (55) 
El{2S - 2Pi/2)exp = 1 057845(9) kHz H, (56) 
El(25- 2Pi/2)cxp = 1057842(12) kHz (57) 



(58) 

Due to large uncertainty and ambiguities with the proton charge radius, one may regard the 
Lamb measurement as a determination of Vp. In this way, from IS Lamb shift, one obtains: 

rp = 0.869(12) fm. (59) 

Logarithmic two-loop corrections significantly alter theoretical predictions for the Lamb shift 
in the single ionized helium as well. The current theoretical value is 

El{2S - 2Pi/2)th = 14 041.57(8) MHz . (60) 



It does not agree with both: the experimental value from |22| and the recent update in |2 
respectively: 

El{2S - 2Pi/2)cxp = 14 042.52(16) MHz , (61) 
EL(2S-2Pi/2)exp = 14 041.13(17) MHz. (62) 

One may wonder about Bqq and further higher order terms, keeping in mind the large value 
of Bqi. There are two possible and complementary undergoing projects: direct calculation 
of this term or numerical calculation of complete two-loop diagrams with Dirac-Coulomb 
propagators. While the second would be the best way, the numerical accuracy might be 
limited at small Z, such as Z = 1. In the direct calculation of Peo one has to consider three 
points: two-loop Bethe logs with e cut-offs, two-loop scattering amplitude with the photon 
mass /i, and the transition terms between e and fi. This project seems to be achievable using 
the methods developed for P50, positronium decay rate and the one apphed here. 
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APPENDIX A: FORMULAS FOR CALCULATIONS OF LAMB SHIFT 

In the calculation of hydrogen and helium Lamb shift we use the following physical constants: 



R = 


10973731.568516(84) 


c = 


299792458 ms'l, 


a'' = 


137.03599958(50) , 


rup _ 


1836.1526675(39), 


me 


rria _ 


7294.299508(16) , 


me 


Tp = 


0.862(12) fm, 


Ta = 


1.673(1) fm. 



(Al) 

In general, Lamb shift in light hydrogen like systems is a sum of nonrecoil, recoil and the 
proton structure contributions. In the nonrecoil limit, known terms are: 
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where is a reduced mass, m = me, and L = \n[m/{fi{Z a^)]. Most of these coefficients 
could be find in any review, such as M or B. The recent result is the direct numerical 



AQo{lS,a) 

Aeoi2S, a) 
Aeo{2Pi/2,a) 

and for He+ {Z = 2) 

Ao(25,2a) 
A6o(2Pi/2,2a) 



-30.29024 
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-0.9735 -0.0640, 
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-30.64466 + 
-0.94940 - 0.0638 





/19 vr^Y 


-0.7961 + 






V45 ~ 277. 



(A4) 
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where the second term is the vacuum polarization p4[. Another recent result is the analytical 
calculation of three-loop contribution in Together with the previously known vacuum 
polarization and anomalous magnetic moment it amounts to 



C40 = 0.417508 . 



(A5) 



In this work we calculate all logarithmic two-loop corrections for S-states. However, for 
P-state only Bq2 is known. For this reason in the theoretical predictions for hydrogen and 
helium we totally neglect higher order two loop corrections for P states, but included B40 
only. We neglect also dependence of in Eq. (pO|) on principal quantum number n, since 
A^ has not yet been calculated for 7^ 1. Recoil corrections, not included in Eq. (|A2|) sum 
to 
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where 



DminSy^) = 4 ln(2) - 



3 - 



2 ' 
/_([+_!) 



(4/2- l)(2/ + 3) 



The finite charge distribution of the nucleus and its self-energy give corrections: 



(A6) 



(A7) 



5E = A(Za)V'r'5w+ ^ 



3n^ 



[Z^a) {ZaY 



In 



M 



/i {Z a) 
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(A8) 



In the theoretical predictions, presented in this paper we have neglected higher order proton 
structure corrections and higher order recoil corrections, which at present are negligible. 
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FIGURES 





FIG. 2. Two-loop diagrams in the Coulomb gauge in NRQED 
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